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Q i Running title: Elliptic Gaudin Model 

C/3 



Abstract. For the elliptic Gaudin model (a degenerate case of XYZ integrable 
spin chain) a separation of variables is constructed in the classical case. The corre- 
sponding separated coordinates are obtained as the poles of a suitably normalized 
Baker- Akhiezer function. The classical results are generalized to the quantum case 
where the kernel of separating integral operator is constructed. The simplest one- 
degree-of-freedom case is studied in detail. 



1 Introduction 



The quantum elliptic (or XYZ) Gaudin model was introduced in [|l|, see also 0, as 
a limiting case of the integrable XYZ spin chain 0. The commuting Hamiltonians 
Hn of the model are expressed as quadratic combinations of s/2 spin operators. 
Determining the spectrum of Hn turned out to be a difficult problem like the original 
XYZ spin chain. Let us list the known facts related to this problem. 

• A solution by means of Algebraic Bethe Ansatz has been obtained only recently 
i. See also §. 

• As shown in [P], in the SU(2)-invariant, or XXX, or rational, case the spectrum 
and the eigenf unctions of the model can be found via an alternative method. 
Separation of Variables, see also the survey [0 . 

• In 1^ the separation of variables in the rational Gaudin model was inter- 
preted geometric Langlands correspondence. 

• In a separation of variables was constructed for the elliptic Gaudin-Calogero 
model which is closely related to the XYZ Gaudin model, though the separa- 
tion of variables for the former one is much simpler. 

• The results of [|] and are based on the interpretation of the corresponding 
Gaudin models as conformal field theoretical models (Wess-Zumino-Witten 
models). The corresponding interpretation of the XYZ Gaudin model was 
obtained in [jlO[, but the conformal field theoretical model corresponding to 
the XYZ Gaudin model turned out to be so complicated that writing down 
the geometric Langlands correspondence for this system, following |p, is not 
easy. 

The main task of the present paper is to present a construction of separated 
variables for the XYZ Gaudin model both in the classical and quantum cases. The 
paper is organized as follows. After giving a detailed description of the XYZ Gaudin 
model in Section ^ we proceed, in Section |^, with the classical case and, following 
the general philosophy of , construct the separated coordinates as the poles of an 
appropriately normalized Baker- Akhiezer function. The corresponding eigenvalues 
of the Lax matrix are then shown to provide the canonically conjugated momenta. 
The whole construction is a simplified version of the one used in P}] . 



The quantum case is considered in Section ^. The separating classical canonical 
transformation is replaced by an integral operator K- We write down a system 
of differential equations for the kernel of K and show that it is integrable. The 
resulting integral operator K intertwines the original and the separated variables 
and provides, respectively, a Radon-Penrose transformation of the corresponding 
D-modules. The quantization constructed is a formal one, since we do not study the 
transformations of the functional spaces of quantum states, leaving it for a further 
study. A detailed study of the spectral problem is given in the simplest case only: 

= 1 (Section |). We show that the corresponding separated equation is none 



1 



other than the (generahzed) Lame equation. Two appendices contain, respectively, 
a hst of properties of elhptic functions, and the formulas describing a realization 
of finite-dimensional representations of s/2(C) on the elliptic curve which are used 
throughout the paper. 

One of the authors (E.S.) is grateful to the Department of Applied Mathemat- 
ics, University of Leeds, UK, where the most part of the paper was written, for 
hospitality and acknowledges the support of EPSRC. The other one (T.T.) is grate- 
ful to the Department of Mathematics, University of California at Berkeley, USA, 
where some part of the work was done, for hospitality. He also acknowledges the 
support of Postdoctoral Fellowship for Research abroad of JSPS. Special thanks are 
to Benjamin Enriquez and Vladimir Rubtsov who gave us their article before 
publication and explained details. 



2 Description of the model 

Let us recall the definition of the XYZ Gaudin model, following The elementary 
Lax operator L{u) of the model depending on a complex parameter u (spectral 
parameter) is given by 

Here a"" are the Pauli matrices, 

. ^ d'udioiu) , . 9[,0oo{u) , . 9[,9o,{u) 

Oil [u) Oil [u) Ooi On {uj 

where Oafiiu) = 6apiu;T), 6^(3 = OapiO), O'n = d/du{9ii{u))\u=o, (see Appendix 0) 
and S"^ are generators of the Lie algebra 5/2 (C): 

[S", S''] = iS\ 

Hereafter (a, 6, c) denotes a cyclic permutation of (1,2,3). Note that A, B, C are 

holomorphic except at u G Z + rZ, where these operators have poles of first order. 

1 2 

Introducing the notation L := L ® ^2 and L := L, where II2 is the unit 
operator in C^, one can establish the commutation relation 

[L (u) , L (v)] = [r{u -v),L{u)+L (v)], (2.3) 
where r{u) is a classical r matrix defined by 

r{u) = -lj2wa{u)cr'^^a''. (2.4) 

^ a=l 
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The r matrix behaves as —■^{V — \)+0{u ^) when -u — 0. Here V is the permutation 
operator: V{x ® y) = y ® x. Exphcitly, in the natural basis in ® C^, 



r{u) 



fa{u) d{u)\ 

h{u) c{u) 

c{u) h{u) 

\d{u) a{u)J 



(2.5) 



where 



a(u) = -b{u) = —, c{u) = , d{u) 

Since Wa{u) are quasiperiodic in u because of ( |A.3| ): 



Wi{u) - W2{u) 



Wi[U) 
W2{u) 
W3(u) 



Wi{u + 1) 
-W2{u + 1) 
-W3(U + 1) 



-Wi{u + r), 
-u?2(m + r), 
W3{u + r), 



(2.6) 



the L operator ( p.l| ) has the following quasiperiodicity: 

L{u + 1) = a^L{u)a\ L{u + t) = a^L{u)a^ 



(2.7) 



Let in {n = 1, • • • , N) be half integers. The total Hilbert space of the model is 
V = 0^^^ Vn, where 14 ^ and 1/^ is a spin i representation space of sl2'- 



p' : shiC) ^ End(c(1^^'^), ^ C2^+^ 

The generating function of the integrals of motion is 



(2.^ 



Tiu) = -TTT\u), 



(2.9) 



where the matrix T(u) is constructed as the sum of elementary Lax operators (|2? 



n=l 



(2.10) 



Here Zn are mutually distinct complex parameters, 



Ln{u) 
oa 



a=l 



1 



Vi 



1 



Vn- 



(2.11) 
(2.12) 



By virtue of the commutation relations (|2.3| ) the operator T(u) satisfies the same 
commutation relations 



[T{u),T{v)] = [r{u~v),T{u)+T{v)], 



(2.13) 
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which imphes the commutativity of f{u): 

[t{u),t{v)] = 0, (2.14) 
Operator t{u) is exphcitly written down as follows: 

N N 

f{u) = pu{u - Zn)L{in + 1) + E ^nCll(« " ^n) + Hq. (2.15) 

n=l n=l 

Here pu, Cii are normalized WeierstraB functions defined by ( |A.5| ) and 

1 3 

rriy^n a=l 
N 3 

Ho = Yl Y.Za{zn-z^)SlS^^ (2.16) 

n,m=l a=l 

are integrals of motion, where 

7 / ,N _ ^10 (^) 7 / ,N _ ^00 (^) 7 / ,N _ ^oi(^) 1 7^ 

^^^'^-i^M^' ^^^'^-i^M^' ^^^'^-i^M^- ^'-'^^ 

Note that the integrals of motion if„ (ra = 0, . . . , A^) appear as coefficients of 



the elliptic Knizhnik-Zamolodchikov equations in |jTO|. Our expression (|2.16| ) for 
Hq differs from that given in because of different normalization of the p and Q 
functions. 

The classical Gaudin model is obtained if we replace all the commutators with 
the Poisson brackets, e.g. 

{T{u),T{v)} = [r{u-v),T{u) + T{v)i (2.18) 

instead of (|2.13|) . The spin variables S"" satisfy, respectively, the Poisson commuta- 
tion relations [S"", S'^] = iS'^ and are subject to the constraint J2a=i{S°')'^ = i"^. 



3 Classical separation of variables 

According to the recipe in [0, the separated coordinates x„ should be constructed 
as the poles of a suitably normalized Baker- Akhiezer function (eigenvector of Lax 
matrix T{u)). The corresponding canonically conjugated variables should appear 
then as the corresponding eigenvalues of T{xn)- Instead of choosing a normalization, 
we shall rather speak of a choice of a gauge transformation M of T{u). The separated 
coordinates Xn will be obtained then as the zeros of the off-diagonal element B{u) 
of the twisted matrix T = M~^TM. 

The classical XYZ Gaudin model is a degenerate case of the classical lattice 
Landau-Lifshits equation for which a separation of variables has been constructed 
TT|, see also a discussion in Here we use essentially the same gauge trans- 



m 



formation M{u) as in |[TT[], and our calculations represent a revised and simplified 
version of those in |TT1 . 
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3.1 Gauge transformation 

Let M{u;u) be the following 2x2 matrix 



M{u-u):=\ 'V ~^'\^yV\, (3.1) 

(^00 [ — ,2) ^00 [ — ,2) 

where u and u are (possibly dynamical) parameters. (This matrix appears also in 
the context of the algebraic Bethe Ansatz. See |T2[, ||^.) A twisted L-operator 
L{u,v] u) depending on a parameter u is defined by 



L(u, v; u) 



^a:^)):-M-.(.;a).,„-„)M,„;.). (3.2) 

Likewise we define the twisted Lax matrix by 

= (e^ifl) -1(;1)) - M-V;a)r(«)M(„;2). (3.3) 

Note that M{u;u) has the quasiperiodicity because of (|A.3| ): 
M(m + 1;m) = -aiM{u;u), 

M{u + t;u) = e-"^("+"/2V3M(M;M)exp(7rm(T3). (3.4) 

These formulae together with ( |2.7|) imply that the function B{u,v;u) has the 
following quasiperiodicity properties: 

B{u+l;u) = B{u;u), B{u + t;u) = e-^^'''B{u;u). (3.5) 



Hence by a standard argument in the theory of elliptic functions (see [T^), we have 



deg(div(5(u))) = 0, -u + J2{multydiv{B{u)))yeZ + rZ, (3.6) 

where multy div(i?(M)) is the multiplicity of a divisor [y] in the divisor dw{B{u)). 
By the definition ( p. 3D , operator B{u;u) is holomorphic except at poles of A{u), 
B{u), C{u), i.e., u = Zn {n = 1,...,N), and at zeros of detM(u;u), i.e., u = 
modulo Z + tZ: 

div(5(u)) > - W + [0]) (mod Z + rZ). (3.7) 
Thus (|3.7|) and ( p.6|) imply that there are (A^ + 1) points Xq, . . . ,xn such that 

N / N \ 

diY{B{u)) = Y.[x,] - N + [0] (mod Z + rZ), (3.8) 

j=0 \n=l J 

and 

N N 

^ = ^ z„ — M (modZ + rZ). (3.9) 

j = 71=1 

Let us fix the parameter u by the condition that one of Xj, for example Xq, is a 
constant ^. Note that u becomes then a dynamical variable. Thus we have 

B{xj;u) = B{u = ^;u) = 0. (3.10) 

Dynamical variables Xi, . . . ,xn are (classically) separated coordinates of the sys- 
tem as we will see below. 
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3.2 Poisson commutation relations and classical separation 
of variables 

The main purpose of this subsection is to prove the following commutation relations. 

Theorem 3.1 Generically the dynamical variables Xj and —A{xj) have the canon- 
ical Poisson brackets: 

(i) {xi, Xj} = for all i, j = 1, . . . , N . 

(ii) {-A{xi),-A{xj)} = for allij = 1,...,N. 

(iii) {—A{xi), Xj} = 6ij for all i, j = 1, . . . , N . 



To prove the theorem, we follow the argument of ||TT|. First let us introduce 
several notations. Define the matrices A, B, C, D as 



-Hj °), J), 0), -I ,3.n, 



Gauge transformation of them are defined as follows: 

A{u;u) := M{u;u)AM{u;u)-^, B{u;u) := M{u;u)BM{u;uy^, 
clu;u) := M{u;u)CM{u;u)-\ D{u;u) := M{u;u)DM{u;u)~\ 



(3.12) 



Bracket (, ) is the standard inner product of the 2x2 matrices: 

{X,Y)=tTXY. (3.13) 

When X{u) is a variable depending on the spectral parameter u, we will denote 
X{xi) by Xi for brevity. For example, 



mem = I 



ti {C{u]u)T{u)). 



U=Xi 



The following statement is proved by the same argument as in the proof of 
Theorem in §2 of [pj]] . 

Lemma 3.2 For any dynamical variable X, 

^^'""^^ ' (3.14) 



and 



We also need the formula for the twisted r matrix. 
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Lemma 3.3 Define 



12 12 

r{u,v;u) := M {u;u)-^M {v;uy^r{u - v)M {u;u)M {v;u), (3.16) 

which we call the twisted r matrix and let rij{u,v;u) be its {i,j) element. Then it 
has the following form. 

fii(M, v; u) = -f22(M, v; u) = -r^isi^u, v; u) = r^^iu, v; u) = 

1 ( e[,{u-v) e[,{u) e[,{v)\ 

2\e,,{u-v) eniu) ^ Oniv) J ' ^ ' ^ 

ri2{u,v;u) = -ris{v,u;u) = -r2i{u,v; -u) = r3i{v,u; -u) = 

= hiiv, u; u) = -r^iu, v; u) = -n2iv, u; -u) = r^^i^u, v; -u) = 
-9'Mv + u) 



29u{u)eniv) ' 
ri4(u, V : u) = r4i{u, v; u) = 0, 

r23{u, V] u) = r32(u, v; -u) - 



9n{u-v)9n{u) 

Proof is given by a direct computation. For example, we have formulae like 



(3.18) 
(3.19) 
(3.20) 



M{u,u)~^(TiM{u,u) 
M{u,uy\ia2)M{u,u) 



1 


^11 


[u)9^i 






1 




^^11 


{u)9n 


[n) 




1 




^^11 


[u)9n 


'u) 



9iq{u)9iq{u) 9io9io{u + u) 

-9io9io{u-u) -9io{u)9io{u) 

9qo{u)9qq{u) 9oo9oo{u + u) 

-9qq9qo{u - u) -9oq{u)9qo{u) 

-9m{u)9oi{u) -9oi9oi{u + u) 

9qi9qi{u-u) 6*01 (?i) 6*01 (m) 



which follow from the addition theorems (cf. [Q, p. 20, p. 22) and the Landen trans- 
formation (cf. |]13| §21.52) of theta functions. Substituting them in the definition of 
f ( p. 16 ) and using the addition theorems again, we can prove the lemma. Q.E.D. 

Proof of Theorem p.l| . 

Using the formulae ( p.l4| ) and ( p.l5| ), we have 
^ idu{CT)),id^{CT))k 



X 



12 12 



id~C T ) 1112 _ _ 

^ " ' ''iC,Ct{T,T}ok) + {C,Ck{T,T},,) 



(3.21) 



12 12 

Therefore computation of {xj,Xk} reduces to computation of {C jC k{T ,T}jk) ■ 
As in Appendix B of [|ll|, we have 



1212 12 1 2 

($,-§,.{T, T}jfc) = tri tr2([$§, f(x,-, Xk)]{f{xj; u) +f{xk] u))), (3.22) 
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for any $, = A,B,C,D. Substituting $ = \1> = C and using {KT^), we have 
X 12 

{CjCk{T,T}jk) = 0. Thus (lOTD imphes that {xj,Xk} = 0. 



A direct consequence of this is {xj, u} = 0, which follows from ( |3.9|) . Using these 
results and Lemma p.2|, we have for j ^ k 



{A{xj),Xk} 



{d^Ckn){A,Co{T,T}^o) - {d~^CoTo){A^Ck{T,T}jk) 
idu{C,T)),{d^CoTo) 



(3.23) 



Hence we need to know {AjC k{T ,T}ji:) and (c^CfcT^). The former can be computed 
by ( ^.22[ ) and ( |3.19| ) and we have 



1112 _ 

{AjCk{T,T}jk) = -2ri2{xj,Xk;u)Ak. 



(3.24) 



The factor (c^CfcT^) is computed as follows: 



{d^CkTk) = {[M{xk;ur%M{xk;u),C]Tixk;u)) 

9[^9ii{xk + u) 



0ii{xk)6ii{u) 



A, 



(3.25) 



Substituting (g:24D, (|3:25|) and ( glSD into (lOl ), we have {A{xj),Xk} = for j ^ k. 

The proof of {A{xj), A{xk)} = is done in a similar way. In addition to the 
formulae we have shown above, we need 



1^12 _ 

{AjAk{T,T}jk) = h3{xj,Xk;u)Cj + ri2ixj,Xk;u)Ck, 

-9[^9n{xk + u) 



29u{xk)ei,{u) 



-Ch. 



(3.26) 
(3.27) 



Proof of the remaining equation {A{xj),Xj} = —1 requires a special care, since 
the r matrix r{u) diverges at m = 0. Instead of (|3.24|) , we use 



1112 _ _ 

{AjC j{T ,T}jj) = -2ri2{xj,Xj; u)Aj - Jim f32(M, Xj; u)B{u; u) 
= -2fi2(xj, Xj] u)Aj + {duB)j. 



(3.28) 



Noting {du{C,T))j = {duB)j and substituting and (lOSi ) into ( |3:23|) , we have 

{A{xj),Xj} = -1. Q.E.D. 

Since B is zero at m = xj, dynamical variable Xj := —A{xj) is an eigenvalue of 

T{x,y. 



fix, 



T{xj 



701 

OOO 1 2 



2 ' 2 



/ 



Xj-\-u r 



701 
000 



2 ' 2 



2 ' 2 



(3.29) 
(3.30) 



Thus if we define the characteristic polynomial by 

W{z, u) := det{z - T{u)), (3.31) 

each pair of dynamical variables {xj,Xj) satisfies an equation 

W{Xj,Xj) = 0, (3.32) 

for j = 1,...,N. Therefore, following the definition in canonical variables 
{xi, . . . ,xn', Xi, . . . , Xjv) are separated variables of the classical ellipitic Gaudin model, 



4 Quantum System: general case 

We return now to the quantum elliptic Gaudin model and construct the quantum 
separation of variables. The special case = 1 is considered in the next section, 
Section ^. 



4.1 Kernel function 

Suppose that the representation space = V^" (|2^ ) is realized as a space of func- 



tions on a certain space with coordinate ?/„ and that the operators 5"^ are differential 
operators on, e.g., polynomials or elliptic functions. The separating operator K is 
expressed as an integral operator 

K f{xi, . . . ,xn) = J dyi- --dyN^ixi, . . .,XN\yi, ■ ■ • , |/iv)/(z/i, • • • ,|/7v), (4.1) 

which maps a function of (?/i, . . . , y^) in Vi (g) ■ ■ ■ (g) Vjv to a function of iV- variables 
Xi on the elliptic curve C/Z + rZ. 

Let us define the operator Xi as follows: 



Lemma 4.1 The following system of partial differential equations satisfies the Fro- 
benius integrability condition. 

B*{xi;u)^ = 0, i = l,...,N, (4.3) 
{Xi + A*{xi;u))^ = 0, ^ = l,...,iV, (4.4) 

where P* is the (formal) adjoint of a differential operator P with respect to {yi, . . . , 
yN) and we set 

N 
n=l 

for a certain constant xq = 



N 



j=0 



(4.5) 
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Proof. This is a consequence of the commutation relation ( |2.13| ). By muhiplying 

1 2 

M{u; u)M{v; u) from the right and its inverse from the left, we have 

12 i 2 

[T{u;u),T{v;u)] = [r{u,v;u),T {u;u) +T{v;u)]. (4.6) 

Note that u is not a dynamical variable in contrast to that in Section |^. 

In order to show the consistency of the equations ( |4.3| ) for i and for j, we 
prove that [B*{xi;u), B*{xj;u)] is expressed as a linear combination of B*{xi;u) 
and B*{xj]u). 

Since the formal adjoint is an algebra anti-isomorphism, (PQ)* = Q*P*, we have 

[B*{xi-u),B*{xj-u)] = [B{xf,u),Bixi-u)]*. (4.7) 

The (l,4)-element of ( |4.6| ) gives 

[B{u] u), B{y- u)] = 2(fi2(M, v; u)B{u) - ¥i2{v, u; u)B{v)) (4.8) 

by virtue of ( |3.19| ) and ( ^.18| ). Replacing u and v in ( ^.8] ) by Xi and Xj respectively 
which are not dynamical, we obtain 

which means that the equation (^4.3|) for i and for j are compatible. 
Next we show the compatibility condition 

[Xi + A*{xi- u),X, + A*{xj; u)] = 0, (4.10) 



which implies the consistency of the equations ( |4.4|) for i and for j {i ^ j). It is 
obvious from (|4.2| ) that 

[X„X,]=0. (4.11) 

Because of ([4. 51), we have 

[Xi,A\xf,^] = ~i^A{x,-u) \ . 



By the same argument as that for ( 3.27 ) the right hand side is rewritten as 



Exchanging i and j, we have 

The (l,l)-element of (|4.6|) means 

[A*{xi]u),A*{xj]u)\ = -ri^{xi,Xj]u)C*{xi]u) - ri2{xi,xfu)C*{xj]u) 

+r3i{xi, Xj] u)B* [xi] u) + r2i{xi,Xj;u)B*{xj]u). (4.14) 
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Summing up ( |4.11| ), ( |4.12| ), ( |4.13| ) and ( [4.14D , we have proved ( [4.10D because of 

dm). 

The consistency of ([4.4| ) for i and (|4.3| ) for j is shown as follows. First assume 
i ^ j- Then the same computation as above gives 



[Xi + A*{xi;u),B*{xj;u) 



'd_ 
du 



B{xj]u)\ +[A*{x,]u),B*{xj;u)] 



Thus we have proved the compatibility of ( |4.4| ) for i and ( [4.3|) for j. Here we used 



(4.16) 



and the (l,2)-element of (f4.6|) . 

The case z = j is almost the same, but there is another term coming from 

[Xi,B*{xi;u)]: 

[Xi + A*{xi;u),B*{xi;u)] = 



d_ 

du 



B*{u;u)- (^B{x,;u)\ + [A*{xi;u), B*{xi;u)]. (4.17) 



By the same computation as (|3.28|) , it follows from the (l,2)-element of ( [4.6[) that 
[A*{xi;u),B*{xi;u)] = 



rB*{xi;u 
On(u) 



d 



nixi)9n{u) ^ du 



B*{u]u). (4.18) 



Substituting (|4l8| ) and (|]T|) for j = i into ( ^TT^ ), we obtain 



[X, + A*(x.;«),S*(x.;S)]= (^-^^^^^^^(x.;^), (4.19) 

\6'ii(u) 6'ii(Xi)/ 

which proves the consistency of (|4.4|) for i and ( |4.3|) for i. Q.E.D. 



4.2 Separating operator 

The separating integral operator K is defined by ( [4.1D with the kernel function 
$(x|?/) satisfying equations ( |4.3|) and (|4.4| ). 

Proposition 4.2 (^zj For any function f of {yi, . . . , z/tv) m Vi ® ■ ■ ■ ® Vat, we /iawe 

ir(5(x,;M)/) = 0, (4.20) 
K{-A{xuu)f) = XJ. (4.21) 

(ii) The elliptic Gaudin Hamiltonian r{u) with the spectral parameter fixed to 
u = Xi is transformed as follows. 



K{f{x{)f){x)=X'^K{f){x) 



(4.22) 



where Xi is defined by U-W - 
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Proof, (i) is a direct consequence of ( |4.3| ) and ( |4.4] ) respectively, 
(ii) By the definition (U) 

if(f(x,)/)(x) = 

^ y <^{x\y) {{2A{xi; uf + B{x^; u)C{xf, u) + C{xi, u)B{xi, u))f{y)) dy 



iA*ix,; u)y^x\y) f{y) dy + J u)B*ix,; uMx\y) f{y) dy 

+ i j [B*{x,- u),C*{x,- u)]^x\y) f{y) dy 



(4.23) 



The first term in the right hand side of ( |4.23 ) is rewritten by the following formula: 

iA*{xf,u)f<l>ix\y) = 



-A*{x,)XMx\y)^ 
Xf^x\y) + [X,,A*{x,)Mx\y), 



where we used ( ^.4| ). The last term of ( |4.24| ) is 



[Xi,A*{xf,u)] 



d_ 

du 



A*{u\u) - 



d_ 

du 



A*iu\u) 



(4.24) 



(4.25) 



because u = ^Zn — Y^Xi. Hence, similarly to the derivation of ( [4.12| ), we can prove 
that 



d 

[X,,A*{x,-u)] = — 



A*(u:u) 



0u0ii{u-Xi) ~ 
2eu{u)9n{xi) 



(4.26) 



The (2, 3)-element of the commutation relation ( [4.6|) gives 
[B {u] u) , C {u] u)] 



2A'(u;u) - .i\ . : B{u;u) 



9[Mu + u) 
Oniu)eniu) 



C{u; u) 



(4.27) 



in the limit v u. Substituting ( |4.24|) , ([4.26|) and (|4.27|) into (|4.23|) and using 
we obtain (^:22|) . Q.E.D. 

The equation ( [4. 201) is a quantum version of ( |3.10D and the equation ([4.21|) to- 
gether with the canonical commutation relation [Xj, Xj\ = 6ij means that operators 
(xi, . . . ,xn', Xi, . . . , Xn) are quantization of the classical separated variables in Sec- 



tion 3.2 



The second statement of Proposition ^]2| provides a formal separation of variables 
for the quantum elliptic Gaudin model. Using the language of and 0, the kernel 
^{x\y) provides a Radon-Penrose transformation of the corresponding D-modules 

(cf. H). 
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In principle, the quantum separation of variables should result in a one dimen- 
sional spectral problem for the separated equation ( [4.22|) which is equivalent to the 



spectral problem for the original Hamiltonians ( p.l6|) . To achieve this goal one needs 
to specify an integration contour in ( ^4.1| ) to study in detail the action of the integral 
operator K on the functional space V . Here we examine only the simplest case 
iV = 1, leaving the general case for the further study. 

5 Quantum system: case = 1. 

In this section we examine the special case of = 1. In this case, everything can 
be computed explicitly and we shall see that the separated equation is nothing but 
the classical Lame equation and its generalization. 

We adopt the realization of the representation of s/2 on the space of elliptic 
functions reviewed in Appendix 0. We could use the standard realization on the 
space of sections of a line bundle over P^, but the result is essentially the same up 
to coordinate transformation and gauge transformation. We omit the suffix n of z„ 
and for brevity. 

5.1 Separated variables 

The quantum twisted B operator B{u] u) = B{u; u) is defined as in the classical case 



(|3J) or (13). Substituting (|2]T|) we obtain 



B{u-u) = fl\ ^(^^10(^-^)^10(^ + ^)5' 

26'ii(m)6'ii(m)6'ii(m - z) ^ 

-9oo{u - z)9oq{u + u)zS^ - 9oi{u - z)9oi{u + u)S^). (5.1) 
The realization of the representation ( p.2|) gives the following expression: 



where 



E(m;m) = 5«(m;u)-^ + 5W(«;m), (5.2) 
dy 



B^'\u;u) = 9u{u)-'9uiu)-'9u{u-z)-'9i,i2y)-' 

/ Z 1i\ f z 

x9,o (^ + --2 + 2) 1^-^+2-2. 



B^'\u;u) = ^^(e,,{u + u)9^,iu-z)9uiyf 

29n{u}9n{u)9n{u - z)9n{yr \ 

+29,0 + + «,„ -y + . - - + -) (-2-2) (5.4) 
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A special point in the case = 1 is that we can make use of the freedom of u 
so that B{u\ u) is a multiphcation operator with the divisor of the form, 

div(5(M; u)) = [x] + [z] - [z] - [0] = [x] - [0] (mod Z + rZ), (5.5) 

as in the classical case, (|3.8|), (|3.9|). In fact, if we put u = —z±2y + l, B^^\u; u) = 
by virtue of (|0| ), and then ( |5.4|) implies 



(We substitute the variable "from the left" , namely we define 

B{u- n)|~_,±2,+i = B^^\u- -z±2y+l)^ + B^'\u; -z±2y+ 1). 

Hereafter we always follow this normal ordering convention.) Therefore we can take 
X = z^2y m (|5.5|) . This is the one of the "separated variables" in this case. 



In the classical model. Theorem |3.1| , —A[x] u) is a dynamical variable canonically 
conjugate to x. This is also the case in the quantum model. The definition of A, 
( |3.3| ), is rewritten in the form 

= Oft I \ft ( AOi,'e,,{u~z)eMe,,{u)s' 

29u{u)9ii{u)9u{u - z) ^ 

-Ooo^ooiu - z)eoo{u)eoo{u)tS' - e,,%i{u - z)eoi{u)eoi{u)s') (5.7) 

by ( |2.1| ). Substituting u = —z±2y + l and u = x = z^2y (from the left), we obtain 

■■z^2y,u=—z±2y+l 

^Ud ( 0M2y) 9'M2y) 9'M2y) \\ 

2\dy ^^10^11(2^) eooen{2y) eoien{2y) ) J 

2\dy p'{y)J ^ ' 

The last equality can be proved by comparing the poles of both hand sides. Therefore 
X = z^2y and X := —A{u] 

'u)\u=z^2yu=—z±2y+i '^^^ the Canonical conjugate variables 

satisfying, 

[X, x] = 1. (5.9) 

We did not make use of the formulation in the previous sections explicitly. In fact, 
thanks to the special choice of u, $ in ( |4.1| ) is a 5-function type kernel, which reduces 
the integral operator K to a. coordinate transformation operator from y to x. 

5.2 Solving the spectral problem 

For the case N = 1, the generating function of the quantum integrals of motion t{u) 
{u is the spectral parameter) 

r{u) = \j:wa{u)\p^'\S'^)f (5.10) 

a=l 
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is explicitly written down. Here we shift the spectral parameter in the original 
definition ( p.9|) as m i— > M2 and set 2; = for the sake of simplicity. Using (p.2| ) 
or ( p.6|) and various identities of elliptic functions in |13|, we can expand the right 
hand side of ( |5.1CI| ): 

(5.11) 



r 



, dy dy^ J 



or 



T 



X 




+ 



2 y?7 — ei 1] 

-l)7] + i{i+l)X 



d_ 
drj 



+ 



(5.12) 



iv - ei){'n - e2){r] - 63) J ' 
where A = p{u). 

As is expected from Proposition O and the result for the rational Gaudin model 
in 0, operator f{u) is factorized as follows when the spectral parameter u is fixed 
to a separated variable xi = 2y. (We may also take xi = —2y.): 



T 



' A. il. 

dy' dy^' 



u 



-^(m;m)L=2j;^= 



-.2y+l) 



(5.13) 



u=2y 



which immediately follows from ( |5.7|) . (The operator X is defined before 
This is consistent with the general result (^4.22|) . 

Equations (|5.11| ) or (|5.12|) shows that the spectral problem of the elliptic Gaudin 
model with = 1 is an ordinary differential equation of second order on the elliptic 
curve C/Z + rZ: 



' d d^ \ 
, dy dy^ J 



(5.14) 



(5.15) 



or on the projective line P^(C): 

^( d d\ 
y' drj' drf' ^ 

Here t(u), t(A) are eigenvalues of f, -i/; G V"'-^^ is an eigenvalue corresponding to 
this eigenvalue. Since operators t{u) and Ua commute with each other by virtue of 
( |B.12| , |B.13| ) and ( |5.10| ), we can decompose each eigenspace of f{u) into those of Ua- 

The equation (|5.14|) has regular singularities: m = (mod F) with exponents 
— 4£, — 2£ + 1, and u = cUq, (a = 1, 2, 3, oji = 1, uj2 = r, uj^ = 1 + r) with exponents 
0, 2£ + 1. The equation (|5.15|) has regular singularities: r] = Ca with exponents 0, 
{2i + l)/2, and r/ = 00 with ^ - -2i. 

If i is an integer, these equations are ordinary Lame equations, while for i G |+Z 
they are generalized Lame equations studied by Brioschi, Halphen and Crawford. 
Following the classical theory of Lame functions (see Chap. XXIII |13|), we can solve 
the spectral problem (|5.14 ), (|5.15|) in V^^^ as follows. 
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5.2.1 Case £ G Z. 

We want a solution ip^r]) of (|5.15|) such that ip{r]) E V^^\ Let us assume that -^(r;) 



is expanded around the singular point as 

oo 

^(^) = E<(^-eaf-^ (5.16) 



r=0 



Oo being 1. The condition ip{r]) G V^^^ means that a" = for r > 2i. Substituting 
(|5.16| ) into ( |5.15|) , we obtain the following recursion relation 



r(£+i-r)a^ = ((£(2£ - 1) - 3(r - 1)(2£ - r + l))e<, + E)a^i 

+ (2£ - r + 2) - r + ^) (e, - efs){ec, - e^)a^_2 (5.17) 

for r > where E = i{i + 1)A — t{X). (Undefined coefficients a'^ for r < are 0.) 
Hence, as a function of E, a" = a"(-E) is a polynomial of degree r of the form 




a^iE) = ArE^ + OiE^-'), A, = j^r! H (^^ - ^ - ^ + J j j (5.1 

Let us denote the roots of a2ij^i{E) = by E'^ {i = 1, . . . , 2£ + 1). The recursion 
relation ( ^.17| ) implies a°(i?f ) = for r > 2£ + 1. Hence we obtain a polynomial 
solution ip{ri) = ip{rj; Ef) of (|5.15| ) of the form (|5.16|) for each i = 1, ... ,2^ + 1, 



provided that 

t{X) = i{i + l)X- E^. (5.19) 

Conversely, ifip{ri) G V^^^ is a solution of the spectral problem ( p.l5| ), then there 
exists certain i for each a = 1,2,3 such that ip{ri) = ip{rj]E^). This is proved by 
expanding the polynomial V'(^) ^is in ( p.l6| ) and tracing back the above argument. 



Proposition 5.1 Assume that uj2 = t is pure imaginary and that parameters Zn 
are all real numbers. Then all Ef are real and the spectral problem ( 5.13^ ) is non- 



degenerate. Namely Ef ^ Ef for distinct i, j and the solutions ip{rj; Ef) span the 
space V^^\ In particular Ef (i = 1, . . . ,2£ + 1) for a = 1, 2, 3 coincide up to order, 
and al^j^-^^{E) = a2£^i{E) = a2£^i{E). Hence we can omit the index a for Ef and 

Vector iplf]] Ei) is an eigenvector of Ua with eigenvalue (—1)^ if af{Ej) ^ and 
(-1)^+1 tfaf{E,) = ^. 

Proof . Under the assumption r G zM, operator f{u) {u G M) is an hermitian 
operator because of ( p.ll|) , and hence it is obvious that Ef are real and that 'ipij]] Ef) 



span V^^\ 

In order to show non-degeneracy of the spectral problem ( |5.15 ) we have only to 
prove that Ef are distinct with each other. Define 



al{E) :- 



al{E), r<£+l, 
'-iy-^al{E), i+1 <r <2i + l. 
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Then the leading coefficients of is 

The recursion relation ( ^.17| ) is rewritten as 

CrCiliE) = qral_i{E) — kr-2^l_2{E) 



(5.20) 



(5.21) 



where 



Cr = r\i+^-r\, 

Qr = (£(2£-l)-3(r-l)(2^-r + l))e, + 
= £-r + ? (2£-r + 2)(ei-e2)(e2-e3). 



(5.22) 



Since ei > 62 > 63 under the assumption of the proposition, we have Cr > and 



kr > 0- This fact together with Ar > (see (|5.20|) ) implies that all the roots of 



al{E) are real and distinct by Sturm's theorem (see, e.g.. Chap. IX, §§4-5, ||T6|). 
This proves the ffist statement of the proposition. 

The operators Ua and f commute and each eigenspace of f is one-dimensional. 
Hence iplrj; Ei) is an eigenvector of Ua- Recall that Ua has eigenvalues (—1)^ with 
multiphcity £ + 1 and (-1)^+^ with multiplicity £. (See §[^3.) If af{E,) ^ 0, then 



-l)V(r/;i?,, 



because of ( P.15|) . Hence there are at most £ + 1 of i^j's such that af{Ei) 7^ 0. In 
other words, at least I of E'j's satisfy a1{Ei) = 0. Since a1{E) is a polynomial of 
degree £, this proves the second statement of the proposition. Q.E.D. 



5.2.2 Case ie^ + Z. 

As in the case £ G Z, we consider an expansion ( |5.16| ) of a solution ip{ri) of the 



spectral problem ( |5.15| ), but this time we consider the series which terminate at 



m = E <iv - e„)''-^ 

r=0 



(5.23) 



They are parametrized by zeros of the polynomial a'^_^_i{E), {-E'f }i=i,...,£+i as in the 
previous case: ip{r]) = tl^^rj; Ef). 

Another set of solutions are obtained from this set by applying the operator Ua'- 



Uai^{mEt) = Y.<iEt){v-^o 

r=0 



(5.24) 



since t/„ and f(u) commute. 



The following proposition is proved in the same manner as Proposition |5Tl . 
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Proposition 5.2 Assume that uj2 = t is pure imaginary and that parameters Zn 
are all real numbers. 

Then all Ef are real and E'f 7^ E'^ for distinct i, j . 

The solutions ^{r]; E"") and Uaip{r]; E") span the space V^^\ In particular E^ 
(i = 1, ...,£+ ^) for a = 1,2,3 coincide up to order, and a],i(E) = a^n.AE) = 

a^_^i{E). Hence we can omit the index a for Ef and 

Vectors ip{ri; Ei) ± Uaip{ri\ Ei) are eigenvectors ofU^ with eigenvalues ^i. 

This proposition means that each eigenvalue E^ degenerates with multiplicity 
two. It was Crawford [|1^ who first found the relation of these two solutions (one is 
obtained from the other by operating U2) by the explicit expansions of type ( |5.23| ), 



([5:2^ ) ■ See also p.578 of fl 



A Notations 



We use the notation for the theta functions with characteristics as follows (see 
for a, 6 = 0, 1, 

0^^(^U]t) = g7rj(n+a/2)V+27ri(n+a/2)(«+6/2)_ ^j^^^ 

Unless otherwise specified, Oabiu) = 9ab{u;T). We also use abbreviations 

d 



9ab = 9ab{0), O'^h - — 



Oab{u). (A.2) 

u=0 



Quasi-periodicity properties of theta functions: 

Oabiu) = i-iyOabiu + 1) = e"™"^,b(M + r). (A.3) 

Parity of thetas: 

Ooo{-u) = 9oo{u), 9oi{-u) = 9oi{u), 9io{-u) = 9io{u), 9ii{-u) = -9n{u) 

A.l Weierstrass functions 

Below we fix cji = 1 and UJ2 = t. 



a{u) = u (1 j exp 

m,n^O ^ 



U I f U ^'^ 



^mn 2 \ CUmn 



(A.4) 



where Umn = + ^^2- 



a'iu] 



a{u) 

(t{u + lji) = -a(n)e'''(2«+'^') 

Ciu + uJi) = Ciu) + 27]i, 

p{u + uji) = p{u). 
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where f]i = C , which satisfy 

Sigma function is expressed by theta functions as follows: 



(t{-u) = -(t{u), C{-u) = -Ciu), pi-z) = p{u). 
Other sigma functions are defined as follows: 



^(^^L+i^) ^00(0) ' 



a(^) ^io(O) ' 



^01 (0) ' 



which satisfy 



Defining ei = p{ui/2), e-i = + a;2)/2), 63 = p(a;2/2), we have 



+ ei = +e2 = + e3 = p{u), 



a^{u) a'^{u) a^{u) 

61 + 62 + 63 = 0. 

ei-e2 = (^)%oi(Or, el-e3=(^)%oo(0)^ 62 - 63 = (^)%io(0)^ 
We also use normalized WeierstraB functions: 

CuH ^ -^Ou{u), Pii{u) ^ --^Ciiiu). (A.5) 



B Realization of spin £ representations on an el- 
liptic curve. 

We recall here the following realization of the spin i representation of the Lie algebra 
s/2(C). Let e, /, h be the Chevalley generators and define = e + f, S"^ = —16 + if 
and 5"^ = h. They satisfy the relation [5", 5*^] = 2iS'^ for any cyclic permutation 
(a, b, c) of (1, 2, 3) and represented by the Pauli matrices cr". 



19 



B.l Spin £ representations. 

The representation space V'^^'^ is realized by 
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k=0 

= { even elliptic function f{y) \ div(/) > -4£(Z + rZ)}. (B.l) 
The generators S"" act on this space as differential operators of first order: 



2 



ru0ui2y)dy e^,iyy 

^ ' e[M2y)dy^ eniyy 

W/C^3^ _ ^01^01 (%) d ^ Myf 

P > e[M2y) dy ^ eniyr 



[B.2) 



or in terms of usual WeierstraB functions, 



(B,3) 



where Ca = p(u;a/2) (a = 1,3,2, uji = 1, U2 = r, = 1 + r) for a = 1,2,3 
respectively and 

_ 1 _ i 1 

Oi — , - , =, 02 — . - . =, 0,3 



(B.4) 

This realization is equivalent to the realization on the space of polynomials of 
degree < 2£ (or, sections of a line bundle on P^(C)), 

e = x''^-2ix, / = -- ^, h = 2x^-2ix, 
dx dx dx 

via a coordinate transformation, x = —9Qi{y]T/2)/9oo{y]T/2), and a gauge trans- 
formation: 



{polynomials in x} 3 <^{x) 



'^)%(.(.)).v.'.. 



Note that this is also obtained by a gauge transformation from a quasi-classical limit 
of the representation of the Sklyanin algebra on theta functions [|1^. 
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The following expression is obtained from the coordinate transformation rj 



l^W = 0Cr7^ (B.5) 

k=0 

and acts on V^^^ as 

pW(5-) = a, (^((e, - e^)(e„ - e,) - {r, - e^f)^ + 2£(r/ - e,)^ , (B.6) 

Let us assume that r is a pure imaginary number. Then, as is well known (see, 
e.g., [0), Ca are real numbers and ei > 62 > 63. This implies that ai and are 
real, while 02 is purely imaginary. 

We introduce the following hermitian form in this representation space: for el- 
liptic functions /(?/), g{y) belonging to V^^'^ defined by (|B.1|) , we define 



{f,9)--= / f{y2)9{yl)^^{yl,y2) (B.7) 

Jc 

where the 2-cycle C is defined by 

C:= {(2/1,2/2) G(C/rf,y2 = yi}, 
and the 2-form /i(?/i,?/2) is defined by 

/^(z/1,2/2) := (ei-e2)^(^+^)(e2-e3)^(^+^) 

a{2y2)a{y2Y'a{2y^)a{y,Y' dy^ A dy, 

''aoo(y2-yi)2(^+^Voo(2/2 + 2/i)2(^+^' 4z 

L ^ (^(2/2) - e2){p{yi) - 62) ^^^^^^^ p'{y2)p'iyi)dy2 A dy^ 

V (ei - e2)(e2 - 63) / 4i ' ' 

This is nothing but a twisted version of the inner product introduced in |T^. If we 
take the description of V^^^ of the form ( |B.5| ), this hermitian form is expressed as 
follows: 



{f,9)--= J^f{r])9{v)Kv,v) (B.9) 
where the 2-form fi{ri, ff) is defined by 

\ (ei - e2)(e2 - 63)7 2i 
An orthogonal basis with respect to this inner product is given by {(?] — e2)-'}j=o,...,2^: 

((r^ - 62^, (r/ - 62)^=) = 27r ^^^||j^^_^~^^^^" (ei - 62)^+^(62 - 63^+^5,,. (B.IO) 

The generators S''' of the Lie algebra s/2 act on the space V^^^-* as self-adjoint opera- 
tors: 

{p'^'\S'^)f,g) = {f,p^'\S'^)g). (B.ll) 

This was first proved in |jl8|, but we can check it directly by using formula ( P.10|) . 

Hence, if u and real numbers, the operator f{u) defined by ( |2.9| ) and 

the integrals of motion Hn defined by (|2.16| ) are hermitian operators on the Hilbert 
space V with respect to (■,■). 
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B.2 Involutions. 

There are involutive automorphisms of the Lie algebra s/2 defined by 

X^{S'') = {-ly-^^^SK (B.12) 
These automorphisms are induced on the spin i representations as 

X^iS') = U-'S'U,, (B.13) 
where operators Ua ■ V^^^ — V^^^ are defined by 

\Vei - e2Vei - 63/ V 2 



for a elliptic function f(y) E V^^\ (cf. [0.) They satisfy commutation relations 

for any cyclic permutation (a, (3, 7) of (1, 2, 3). The action of these operators on the 
bases {(// - e„)^}j=o,...,2£ is: 

U2iv-e2y = e'''^''-^\e,-e2y-\e2-e,y-'ir]-e2f-\ (B.15) 
Usiv-esY = e--^'{e,-esy-'{e2-e,y-\r]-esf-^. 

Hence eigenvalues of Ua are (—1)^ with multiplicity i + 1 and (—1)^+^ with multi- 
plicity £ if £ is an integer, and ±i both with multiplicity i + ^ if ^ is a half of an odd 
integer. 

When LJi = 1 and 002 is a pure imaginary number, these operators are unitary 
with respect to the hermitian form 
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